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Recently, we have experimentally demonstrated a continuous loading mechanism for an optical 
dipole trap from a guided atomic beam [1]. The observed evolution of the number of atoms and 
temperature in the trap are consequences of the unusual trap geometry. In the present paper, 
we develop a model based on a set of rate equations to describe the loading dynamics of such a 
mechanism. We consider the collision statistics in the non-uniform trap potential that leads to two- 
dimensional evaporation. The comparison between the resulting computations and experimental 
data allows to identify the dominant loss process and suggests ways to enhance the achievable 
steady-state atom number. Concerning subsequent evaporative cooling, we find that the possibility 
of controlling axial and radial confinement independently allows faster evaporation ramps compared 
to single beam optical dipole traps. 
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I. INTRODUCTION 

For almost two decades, the standard way of cooling 
atoms to ultra-cold temperatures has been: laser cooling 
in magneto-optical traps followed by evaporative cooling 
in magnetic and optical traps. For many atomic species 
and especially molecules, different methods are required 
since laser cooling on a closed cycling transition is tech- 
nically not feasible [51 [3]. In recent years, there has been 
considerable progress in the generation of slow and inter- 
nally cold atomic and molecular beams [IHS], in particu- 
lar for atoms or molecules with an electric [THIO] or mag- 
netic [TTJ [T2] dipole moment. A mechanism that allows 
for a further increase of the beam's density p!5H2D] may 
hence provide suitable starting conditions for evaporative 
[21j or demagnetization cooling [22 . Such a mechanism 
can be realized in a purely statistical way [17' or by scat- 
tering a single photon 23J that changes the internal state 
from untrapped to trapped. 

We have recently demonstrated a loading scheme that 
is based on the accumulation of chromium atoms from a 
continuous flux of guided ultra-cold atoms pH - f2S| in a 
conservative potential and allows for the fast production 
of Bose-Einstein condensates (BECs) [T]. This scheme 
benefits from the strong magnetic moment of chromium 
atoms. It is likely to work also for other species with a 
comparably versatile level structure and strong magnetic 
moments such as atomic erbium |27| or dysprosium |28j 
that have recently gained interest. The scheme might 
even be the favorable technique of producing BECs from 
such elements since it omits some of the limitations that 
are generic to those atoms, e.g. strong dipolar relax- 
ation losses and light assisted collisions. In the case of 
chromium, the repetition rate has increased from an ex- 
perimental cycle of more than a minute [29 to 5 seconds 
for the production of a BEG [1]. 

In earlier work [55] , we estimated the transfer efficiency 
from the guided beam into the trap from the evaluation 
of single particle trajectories. That approach allowed a 
prediction of the loading rate and agrees with our mea- 



surements. However, we observe a saturation of the atom 
number at an unexpectedly low level. The estimated 
maximum number of atoms due to three body collisions 
of ground state atoms or light assisted collisions should 
be an order of magnitude higher than the observed ones. 
In this article, we develop a model based on coupled rate 
equations to quantitatively describe the loading mecha- 
nism. With our model we are able to identify elastic two 
body collisions, i.e. evaporation as the limiting mech- 
anism in our current trap. The model takes the trap 
potential's specific non-uniform geometry into account. 
Regarding evaporation, this implies a dimensionality [30] 
of two. Unlike one and three dimensional geometries [31] . 
this has to our knowledge not been treated so far. 

The paper is organized as follows: In section II we first 
describe briefly the experimental procedure used and dis- 
cussed in detail in jH [24H26] . We develop a rate equation 
model for the loading process in section III, and com- 
pare numerical simulations based on the model to our 
experimental data in section IV. 



II. EXPERIMENTAL SETUP 

In the following, we briefly summarize the underlying 
mechanism of the loading scheme for convenience [U I20j . 
Fig. [l] shows a schematic illustration of the experimen- 
tal setup. We overlap a single beam optical dipole trap 
(ODT) on axis with a continuous beam of ^^Gr atoms. 
Thereby, about half the flux of atoms in the beam is 
radially confined in the optical potential. In order to dis- 
sipate the atoms' directed kinetic energy, we place a pair 
of coils around the focus of the ODT to create a magnetic 
field maximum along the atoms direction of propagation. 
Since the atoms are initially in a low-field seeking sub- 
state, they are slowed as they approach the magnetic 
barrier potential. Glose to their classical turning point, 
an optical pumping beam transfers the atoms to their ab- 
solute ground state, i.e. a high-field seeking state which 
is attracted axially by the magnetic saddle potential ere- 
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FIG. 1. (Color online) Schematic illustration of the setup un- 
derlying the continuous loading mechanism. Atoms are trav- 
eling to the right and are funneled into the radially confining 
ODT potential. Axially, the barrier coils create a state de- 
pendent potential that is switched from a wall to a well by 
means of the optical pumping beam. 

ated by the barrier coils and radially confined by the 
(ODT). We use the following experimental parameters: 
the ODT is focused to a waist of 30 /im and operated 
at a power of 80 W which results in a potential depth of 
about 1 mK. Its focus coincides with the geometric center 
of the magnetic barrier coils. They have a diameter of 
1 mm, are spaced at a distance of about 1 mm and are 
driven by a current of 1 A. The atomic beam has a ve- 
locity of Im/s. All experimental information presented 
in this article is obtained from in-trap and time of flight 
absorption images taken along the vertical axis. 

III. MODELING THE LOADING PROCESS 
WITH RATE EQUATIONS 

In the following section, we develop a model based on 
rate equations to explain the physics of the loading pro- 
cess. We will later use the model to interpret our experi- 
mental data. Starting from an empty trap, the evolution 
of the atom number is proportional to the loading rate 
L until loss processes that depend on the atom number 
become dominant and eventually an equilibrium state is 
reached. These considerations, lead to the common rate 
equation for the evolution of the atom number N: 

N^L-^N- -fevN^ (1) 

where 7 is the single-body loss rate that includes e.g. col- 
lisions with molecules from the background gas and fast 
impinging atoms from the guided beam. accounts for 
losses due to binary collisions between trapped atoms. 
Since the atoms are in their electronic ground state and 
dipolar relaxation is forbidden due to energy conserva- 
tion [32j , only elastic two-body collisions can occur in our 
system. During elastic collisions selectively hot atoms es- 
cape which leads to a net cooling of the trapped sample. 
On the one hand, the last term in Eq. ([!]) thus also 
changes the temperature, while on the other hand j^v it- 
self strongly depends on the temperature of the trapped 



gas. The model described by Eq. ([T]) hence needs to 
be extended to a set of differential equations that also 
takes into account the evolution of the temperature. In 
the following section we deduce this model from basic 
considerations 



A. 2D-Evaporation 

The physical picture behind evaporation is based on 
the distinction between trapped and untrapped momen- 
tum states |31j . A particle initially occupying a trapped 
state can be promoted to an untrapped one through colli- 
sions with other particles. We restrict our considerations 
to the Knudsen regime (mean free path larger than the 
size of the cloud), i.e. every particle promoted to an un- 
trapped state escapes. 

The dimensionality of evaporation is given by the 
condition of truncation in momentum space |30| . 3D- 
evaporation is the simplest case: a particle escapes un- 
der the condition e > ej, i.e. the particle's total energy 
e = U {r)+p'^/ {2m) exceeds the trap depth e^, where U{r) 
is the potential energy at the position r and /{2m) the 
particle's kinetic energy. 3D-evaporation occurs in po- 
tentials of uniform depth in all directions (e.g. ideally a 
single beam ODT or a magnetic trap with rf-knife). 

The loading mechanism discussed here is based on the 
superposition of an ODT with a magnetic field maximum 
[D i2QJ . Geometric constraints imply that the hybrid 
trap potential is to a good approximation separable into 
an axial (purely magnetic confinement) and a radial part 
(purely optical confinement): 

U{p,z) = U{p) + U{z), (2) 

where the potential depth in the radial direction Cp is 
much smaller than along the axial direction et.z for typi- 
cal experimental parameters. Therefore, the criterion of 
evaporation is more strict: the particle's energy in the 
radial degrees of freedom must exceed the radial trap 
depth Up{p) + Pp > €t,p- This reduces the probability 
for hot atoms to be removed and thus evaporation is 
less efficient. 3D-evaporation, however, can also occur 
in a non-uniform potential trap due to ergodicity |37) . 
Whether or not evaporation can be considered ergodic 
depends on whether a particle finds the trap's exit suf- 
ficiently fast in the time between successive collisions. 
For the trap potential discussed here, numerical evalua- 
tion of the classical single particle trajectories shows that 
the dynamical coupling between the degrees of freedom 
is negligible. Surkov [5D] and Pinkse [SI] compare ID- 
evaporation and 3D-evaporation based on an experiment 
on magnetically trapped atomic hydrogen. For the de- 
scription of the loading process discussed here, a model 
considering 2D-evaporation is required. 



3 





ative momentum q = {pi — p2)/2. Energy conserva- 
tion requires the relative momentum's absolute value 

q = \q\ ^ \f<il 



FIG. 2. (Color online) (left) The intersection of the surface S 
of a sphere with radius q and a long cylinder with radius Qp 
geometrically illustrates the loss criterion for particles. The 
distance between the center of the sphere and the cylinders's 
axis represents Pp. The fraction of S that lies outside the 
cylinder corresponds to the statistical probability /loss of one 
(or both) atoms to escape the trap, (right) cut along the px 
Py-plain. Example: The q vector, pointing along the p^-axis 
before collision, has reoriented into the p^ Py-plain and both 
particles now occupy an untrapped state. Alternatively, only 
one atom would be promoted to an untrapped state, if the 
momenta was reoriented along the p^j-axis. 



B. Collision statistics 

In this section, the Ansatz for ID-evaporation treated 
in |31| is extended to the case of a non-uniform poten- 
tial with cylindrical symmetry. The most convenient way 
to describe the elastic collision of two particles with mo- 
menta pi and p2 is to use center of mass coordinates 
with the average momentum p = (pi + P2)/2 and rel- 



qy + q1 to be conserved. Assuming 

s-wave collisions, a collision thus reorients q isotropically 
on the surface of a sphere with radius q. Furthermore, 
each component of the average momentum p is conserved 
during collisions. With the assumption et.z ^ Ct.p and 
the cylindrical symmetry of the trap potential (Eq. ([2])), 
particles leave the trap if their radial momenta pp after 
the collision are larger than the radial escape momentum 

Qp- 

Within the center of mass coordinates, it is merely influ- 
enced by two parameters whether a collision redistributes 
particles into trapped or untrapped states: q and Pp. q 
denotes the relative momentum of the colliding particles, 
which is randomly reoriented (azimuthal and polar angles 

(/>, 9) after a collision and Pp 



Ip^. + P^, the average ra- 
dial momentum, respectively. Geometrically, the escape 
condition thus equals the intersection of a sphere's sur- 
face (radius q) with an infinitely long cylinder (radius 
Qp) with an offset of Pp as depicted in Fig. [2| States 
outside the cylinder are untrapped. It should be noted 
that, in contrast to the case of 3D evaporation, an elastic 
binary collision in a non-uniform potential can also lead 
to the loss of both particles. 

Mathematically, Eqn. ([3])-(|4| express the probability 
/loss of one (or both) atoms to escape the trap as a 
function of the input momenta q and Pp by integration 
over all possible reorientations of q. 
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pI) sin 9 d9: 



(3) 
(4) 
(5) 
(6) 



where B(pp,q,9,4)) is the criterion for an atom to leave 
the trap. & represents the Heaviside step function. The 
mean energy per particle in the radial ejjj" = + and 
axial e+ degrees of freedom after the collision is obtained 
by integration over all possible reorientations of momenta 
as stated by Eqn. ([5])-(|6|. 



C. Thermodynamics 

Eqns. (|3])-(|6]) express the statistical results for general 
collisions in the trap as functions of the input momenta 
q and Pp. The goal, however, is to write down /loss, , 
and e+ depending on temperature. This is achieved by 
summing Eqn. over the statistical distribution of 

q and Pp at a given temperature. For a thermodynamic 
description of the collision statistics we make the follow- 
ing assumptions: 
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• The axial and radial degrees of freedom individu- 
ally are in quasi equilibrium while different temper- 
atures along the respective axes are allowed. 

• The thermal momentum distribution f{pi.z,Tz) of 
a single particle along the axial direction is given 
by a Boltzmann distribution. 

• The thermal momentum distribution f(pi^p,Tp) is 
given by a truncated Boltzmann distribution |35| . 

• On time average, equipartition of mean kinetic and 
mean potential energy is assumed (as for the case 
of a harmonic potential). 

The first three points are a direct result of the trap's 
prolate symmetry (regarding trapping frequencies and 
depth) in combination with the atomic beam. The atoms 
loaded into the trap have different temperatures in ax- 
ial and radial directions. These assumptions now allow 
to describe the atomic cloud with two temperature pa- 
rameters: Tz and Tp. Experimentally, both values can 
be accessed by time of flight images. The momentum 
distribution of a single particle within an ideal gas at 
equilibrium is given by the Boltzmann distribution: 



f {Px,Py,Pz 



1 



2nm,kBT 



3/2 



(7) 



A gas trapped in a potential with a finite depth can be 
approximated by a Boltzmann distribution truncated at 
the trap depth e^. It is then convenient to introduce the 
truncation parameter 77 = et/fcsTp = Q^/2mkBTp. To 
derive the statistical distributions for q and Pp at a given 
temperature, Eq. ([t]) is transformed into the respective 
coordinate systems, which leads to (for a derivation see 

my- 

fq{q,Tp,Tz) 



no Pp e-fp/™'=«'^'' 



(8) 



riQ gerf q. 



Tp~Tz 
mkBTzTp 



-q ImkuTp 



(9) 



Fig. |3] shows examples of these two thermodynamic dis- 
tributions. In Fig. |3](a), fpp{Pp,Tp) is plotted for differ- 
ent trap parameters 77. The distribution is truncated at 
the escape momentum Qp. fq{q,Tp,Tz) is drawn in Fig. 
[3] (b). The two curves compare a cross-dimensionally 
thermalized distribution (black) with a non-equilibrium 
one (gray) at an equal mean energy. We sum over all 
possible momentum states by integrating Eq. |4]) multi- 
plied by their statistical weights Eqn. ([8| and Q, respec- 
tively. This yields the desired function for the probability 
of evaporation [55] : 

/loss (^p ; Tz ) — 

= n-t / fppiPp) fqiq) fiossiPp,q) dppdq (10) 

JO Jo 




FIG. 3. (a) Examples of the statistical distribution of mo- 
menta fpp{Pp,Tp)\^^^Q2/2mkBV © different 
values of rj. (b) Plots of the statistical distribution of the mo- 
menta fq{q,Tp,T2) given by Eq. (|9|. The black and the gray 
curves describe particles with the same mean energy e. The 
black curve shows the thermodynamic equilibrium {Tp — Tz), 
while the gray one represents the case where Tp <Tz. 



with the re-normalizing pre-factor nt{Tp) = 

ilo^" fppiPp) dpp)^^. In the same way the distri- 
butions of the kinetic energies (per particle) among the 
axial and radial degrees of freedom after a collision are 
obtained: 



00 i-Qp 

nt I j fppiPp) fqiq) epiPp,q)dppdq 

(11) 

^{Tp, T^) =nt I I fp^ (pz) fq{q) et{pz,q) dpz dq, 
Jo Jo 

(12) 

where the integration along pz simply adds a term of 
(2 — fioss)kBTp/2. Another consequence of the finite ra- 
dial trap depth is that ip ^ kBTp, i.e. the mean ther- 
mal energy Ip per particle is no longer proportional to 
the temperature parameter Tp. This can be easily seen 
from Fig. [sjja). The mean energy is thus a function of 
the temperature parameter and the escape momentum. 
Conversion between ep and Tp is done with the following 
formula and its numerically evaluated inverse: 



(-p{Tp,Qp) = nt{Tp) 



Qp 

2m 



fppiTp, Pp) dpp (13) 



The obtained f\ass{Tp,Tz) is plotted in Fig. |4] One can 
see that (for high temperature parameters) on average 
more than one particle can be lost after a binary colli- 
sion. The averaged loss or gain of kinetic energy after one 
collision can now be written as the difference between the 
mean thermal energy before e and after a collision : 



(14) 
(15) 



The considerations so far and in particular Eqns. (14) 



and ( 15 1 only include kinetic energy, i.e. they apply to 
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FIG. 4. The function fioss[Tp,T^) defined in Eq. ([To| is plot- 
ted with respect to generalized coordinates. The x-axis cor- 
responds to the truncation parameter. Similarly, the y-axis 
relates the radial trap depth to the axial temperature. 



a trap shape where potential energy can be neglected, 
as in a box of finite size and infinitely steep walls. In 
a harmonic potential, however, a particle on time av- 
erage has an equal amount of kinetic and potential en- 
ergy Cpot = Ekin- During an impact, the kinetic energy 
Ekin = i^/2m is redistributed among the colliding parti- 
cles as described by the model derived above, while both 
colliding atoms retain their fraction of potential energy 
before the collision. To account for the potential shape, 
we assume in the following that only half of a particle's 
total energy can be altered in a single collision event: 



^^P "n'^p + 9 (2 - /loss) £p ~ Cp - ^ ~ /loss Cp) 



(16) 



=\4 + \{'^- /loss) f-z - 



/loss 



(17) 



The shape of the potential determines the distribution 
between kinetic and potential energy [3D] and in this way 
it is responsible for the speed of cross-dimensional ther- 
malization. 



D. Collision rate 



which scales with the product of average density, scat- 
tering cross section cr, and the mean relative velocity. 

To obtain the trap's effective volume for evaporation, 
the phase space distribution of a classical ideal gas at 
thermal equilibrium j35^ 



/(r,p) = 710 



2Trh 
TTmkBT 



\U(r)+p^/2m\/kBT 



(18) 



is integrated over all momentum states that occur in the 
trap potential: n(r) = (27r?i)~^ / cPpf {r,p). With the 
discussed properties of the potential (separability and 
non-uniformity) this leads to: 

n(r) = no (l - el^''-^-^"'!/'^-^-) e^t^.A-^p • e'^^/'^-^^. 

(19) 

The effective trap volume of the trapped gas for evap- 
oration is defined by ^33) : 



J 71^ (f) (Pr 



(20) 



For cr, we use the energy dependent s-wave scattering 
cross section in the effective range approximation |42) : 



a{k) 



fc2a2 -I- (l/2fc2are-l) 



(21) 



with a = 103ao and the effective range Tg — 74ao [43] . oo 
denotes the Bohr radius. Contributions of higher order 
partial waves are neglected although they are likely to 
play a role at the present temperature range. However, 
related experimental or theoretical data are lacking |44j . 

We define the collision parameter 7coii through the re- 
lation [551: 



Tcoii = Tcoii N 



(22) 



This leads to the definition of the thermally averaged 
collision parameter: 



j;^2qa{2q)fgiq,Tp,T,)dq 



(23) 



The relative wave vector k in the scattering cross section 
is substituted by the relative momentum coordinate q — 
hk, where the factor of 2 in front of q stems from the 
relative mass /z = m/2. 



In Sec. nil CI we derived the statistical outcome of one 
single collision. In the following we assume that the evo- 
lution of the trap (considering N, Ep, and E^) is deter- 
mined by the product of the collision rate and the mean 
effect of a single collision [41]. The rate at which these 
events occur is given by the elastic collision rate Fcoii 



E. The rate equations 

Together, the considerations from the previous para- 
graphs result in the following set of coupled differential 
equations: 
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N ^ L ~ 7con(T;,r,)/ioss(T;,T,)iV2 

Ep = EL, pi " l^pN - 7cou(Tp,T'z) • [e^{Tr,T^) 



(24) 

fioss{Tr,T,)ep]N^ (25) 
/ioss(rr,T,)e,] (26) 



with 



2e. 



2E, 
Nkf 



Tp^Tp{ep)^Tp{^ 



(27) 
(28) 



where Eq. (28) is the inverse function of Eq. (13). E 



and Ez are the total kinetic energies in the radial and 
axial degrees of freedom, respectively. The first terms on 
the right-hand side of Eqn. ( 24 )-( 26 ) describe the loading 



process. Particles are inserted into the trap at a loading 
rate L, a mean kinetic energy of Il,z axially and eL,p 
radially. Single body loss processes such as background 
gas collisions and fast atoms from the beam are covered 
by the second terms with the respective coefficient 7. 
Finally, the last terms account for binary elastic collisions 
among trapped atoms. They include plain evaporation 
and cross-dimensional thermalization within the model 
developed in this section. 



IV. THE EXPERIMENT 

A. The loading process 

Fig. [sjja) shows the evolution of the atom number, 
temperature, and PSD during the loading process. At 
time t — 0, the optical pumping beam is switched on and 
thereby the loading process starts. Experimental data 
extracted from in-trap absorption images (squares, trian- 
gles, and circles) are compared to the model developed 
in Sec. Ill The theoretical curves (solid lines) are based 
on the numerical integration of Eqn. (|24|) 



The set of 

coupled differential equations require three starting val- 
ues: iVo, Ep Q, and E^^- Since the loading process starts 
from an empty trap, these parameters are all zero. Fur- 
thermore, the equations contain the parameters 7, eL,p, 
eL,z, and L. 7 is estimated from decay measurements 
of the fully loaded trap in presence of the atomic beam. 
It lies well above 1 s and therefore does not play a role 
for the loading process under typical experimental con- 
ditions. The remaining parameters are directly deduced 
from the respective experimental data, within their re- 
spective experimental uncertainty. L is obtained from a 
slope of iV(i) in the initial stage. £l,p and eL,p correspond 
to the temperature parameters Tp and (within their 
experimental uncertainties) of the cloud before collisions 
play a role. 



One can see that during the first 50 ms of loading, the 
atom number increases at a constant rate. A linear fit 
yields a loading rate of 5.8 • lO^atoms/s, which corre- 
sponds to a loading efficiency of 10-15% with respect 
to the flux in the atomic beam. After about 100 ms N 
saturates around N = 5 ■ 10^ atoms. As also shown in 
Fig. [5](a), the temperature decreases during the loading 
process. 



B. Plain evaporation 

As a benchmark of the loading scheme it is interesting 
to validate the agreement between the model and exper- 
imental data under different conditions, i.e. after the 
loading process is terminated and thus L — 0. In Eqn. 
(24)- (26) the first terms on the right hand side vanish. 
7 is smaller than during the loading process. Since the 
atomic beam is switched off, the 1 /e lifetime is now solely 
limited by background pressure to 30 s. The remaining 



starting conditions for the simulation are: Nq, Ep Q, and 
Ez.o- These are obtained from measurements that are 
taken directly after the loading process has been termi- 
nated. 

Figs, [sjb) and (c) show the temporal evolution of the 
trap directly after the loading process has been stopped 
with a preceding ijoad =100 ms (b) and tioad s (c), re- 
spectively. Data collected during plain evaporation are 
compared to the model developed in Sec. Ill The plots 



in the top row show the evolution of the atom number, 
while below the respective temperature measurement is 
displayed. The plots cover the first two seconds, which 
is the experimentally most relevant time frame. To em- 
phasize the dynamics at small and large time scales, the 
insets have an extended x-axis up to 30 s and a double 
logarithmic scale. 

From the evolution of N in Figs. [5f^b) and (c), one can 
see a fast decay during the first 500 ms that can be at- 
tributed to plain evaporation before eventually, the decay 
becomes exponential. As expected, the initial fast de- 
cay is more pronounced in (b), where the trapped atoms 
have a higher mean energy compared to (c). Concerning 
temperature, one observes Tp and Tz equilibrate within 
about 50 ms, as well as a decreasing mean temperature 
throughout plain evaporation. Despite the atom losses, 
the PSD is increasing by about one order of magnitude 
during plain evaporation, which indicates that evapora- 
tive cooling due to elastic collisions is the dominant pro- 
cess. 
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(b) Plain evaporation after t,„„j= 100 ms 
xio' 
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(c) Plain evaporation after t,„,j= Is 
xlO^ 
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FIG. 5. (Color online) Evolution of the number of atoms A'^, the temperature parameters (Tp and Tz) and the PSD during 
the loading process (a) and successive plain evaporation after fioad = 100 ms (b) and tioad = 1 s (c). Data are extracted 
from in-trap absorption images (squares, triangles, and circles), while the solid lines show computations according to the model 
developed in Sec. |III| Displayed errorbars result from an estimated uncertainty in the magnification of the imaging system of 
10%. 



C. Limitations to the loading process 

As can be seen from Fig. [5j the number of atoms 
saturates after a trap loading time of less than 200 ms. 
It would be favorable to prolong the loading phase and 
reach a higher equilibrium atom number iV^ax- There- 
fore, in the following, we investigate the mechanism that 



is responsible for the saturation of N. The good agree- 
ment between experimental data and simulations sug- 
gests plain evaporation to be the dominant loss process. 
No other losses are included into the model, apart from 
single body losses, which are yet on another time scale. 
Fig. [6]Ja) shows simulations and experimental data of the 
loading process with different loading rates. We exper- 
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ior described by Eq.(29). It might seem counter-intuitive 
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FIG. 6. (Color online) (a) The temporal evolution of the 
atom number is plotted for data corresponding to different 
loading rates L. The squares represent experimental data 
and are compared to simulations (solid line), which merely 
differ in the loading rate L that is measured directly from 
the initial slope of the experimental data, (b) The squares 
show the equilibrium atom number A'^max measured after a 
long fioad with respect to the corresponding loading rate L. 
The diamonds connected by a solid line show A'^max from the 
simulations. 



imentally reduce L by decreasing the flux of Cr atoms 
between our effusion cell and the MOT. The measured 
loading rate is then used in the respective simulations, 
while all other parameters remain the same. When an 
elastic two-body loss process such as evaporation consti- 
tutes the dominant source of loss, the stationary solution 
of the simplified rate equation for the atom number Eq. 
([T]) is of the following form (i.e. ^evN ^ 7): 



N — 




(29) 



is plotted as a function of L in Fig. |4][b). Both 
simulations and experimental data agree with the behav- 



that elastic collisions between trapped atoms can actu- 
ally limit A'fnax- Atoms are loaded into the trap with a 
mean energy corresponding to 77 « 2. The mechanism can 
therefore be understood in the context of the continuous 
heating rate imposed by the loaded atoms and the cool- 
ing rate due to evaporation. An equilibrium temperature 
develops where jco\i{Tr,Tz) ■ fioss{Tr,T^) = 72 = const. 
The key to increase -/Vmax, is therefore a reduction of /loss- 
One way to achieve this is to increase the radial trap 
depth, e.g. by a higher ODT power. Simultaneously in- 
creasing the ODT's waist maintains the trap frequencies 
and thus avoids additional radial heating. Our compu- 
tations predict that loading the trap at 77 = 3 results in 
gain by a factor of 2 in iV^ax compared to the current 
experimental parameters (at equal L). At the same time, 
the steady state PSD improves by more than a factor of 
7. 



D. Forced evaporation 

Plain evaporation as described by the model developed 
in Sec. |III| quickly slows down with increasing rj. To 
maintain evaporative cooling, the trap depth has thus 
to be continuously lowered (forced evaporation). In our 
case, we realize this by decreasing the ODT's power. 
There is a large range of literature on evaporative cooling 
in magnetic traps l^niEIlEllllS] or in optical dipole traps 
[IB, ,47^ . The hybrid magnetic and optical trap geome- 
try, as a prerequisite of the continuous loading process, 
combines advantages of either trap type. On the one 
hand, the high-field seeking absolute ground state can 
be trapped, while on the other hand, the trap weakens 
the typical problem of a single beam ODT: the decreas- 
ing collision rate as a result of the coupling between the 
trapping frequencies and trap depth. Here, the axial con- 
finement is maintained throughout the evaporation ramp 
by the magnetic field, while merely the radial trapping 
frequency is lowered. In order to extend the model de- 
scribed above for forced evaporation, we have to consider 
the work exerted on the trapped g result of the 

variation of the radial trap shape |31| : 



NksT f dVe \ dPoDT 



K(Tp,T,) \dPoDT 



dt 



(30) 



Spilling dS], i.e. the loss of trapped states due to changes 
in the potential shape can be neglected for our experi- 
mental conditions (77 7 during the ramp). 

Fig. [T^a) shows the experimentally optimized ramp of 
the ODT power for forced evaporative cooling. It consists 
of one plain evaporation segment starting at t = and 
subsequently 5 segments with linearly decreasing ODT 
power. In (b), the evolution of the PSD is shown with 
respect to N. The squares show data taken during the 
evaporation ramp, while the solid line indicates the re- 
spective prediction according to model of the rate equa- 
tions. It should be noted, however, that these do not 
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FIG. 7. (Color online) (a) Experimentally determined ramp 
of the ODT power that is used to produce a BEC. At t = the 
loading process is stopped and followed by plain and forced 
evaporation as indicated, (b) Evolution of the PSD as a func- 
tion of the atom number during evaporation. The squares are 
data points that are measured during the evaporation ramp. 
The error bars result from fit uncertainties and systematic 
errors as the magnification of the imaging system. The solid 
line represents a simulation of the forced evaporation process 
according to the classical model discussed in the text. 



consider quantum statistical effects. After forced evapo- 
rative cooling, we reach a BEC with N « 10'' atoms in 
an overall time of 4.5 s. 



modeling continuous loading schemes for an ODT from a 
beam of atoms [U [20] . Due to the non- uniform trap ge- 
ometry and the simultaneous appearance of loading and 
evaporation, the dynamics of the thermodynamical prop- 
erties of the trapped cloud are far more complicated than 
in previously modeled systems. As a direct result of the 
implementation of the loading scheme, the dimensional- 
ity of evaporation [30] is two for the non-uniform trap 
geometry discussed here. 

We have shown that numerical simulations of the load- 
ing process, plain, and forced evaporation agree quan- 
titatively well with experimental data. In experiments 
as well as in numerical simulations we observe the onset 
of saturation at an atom number of about 5 • lO"" after 
loading the trap for 100 ms. The model allows identify- 
ing elastic binary collisions as the limiting process within 
our current experimental conditions. It is important to 
remark that the system is not yet limited by light assisted 
collisions mediated by (re)absorbed pump-light photons. 
Larger steady state densities can therefore be reached 
through a higher trap parameter rj of the atoms loaded 
into the trap. This may be realized by means of a radi- 
ally colder beam or a higher radial trap depth. 
Furthermore, we have shown that the specific trap geom- 
etry has good properties for successive evaporative cool- 
ing, since the axial trap frequency is independent from 
the trap depth. In combination, the loading mechanism 
followed by forced evaporation in the hybrid magnetic 
and optical trap allows producing a Cr BEC within less 
than 5 s, which is significantly faster than previously pub- 
lished [151 HH] experiments on Cr atoms. 
The demonstrated loading mechanism may as well be 
applicable to atoms where a closed cycling transition for 
laser cooling is not available, or for molecules, which can 
except for a few cases not be laser cooled due to their 
complex level structure. Thereby, good starting condi- 
tions for successive evaporative cooling can be reached, 
circumventing the need for laser cooling on a closed cy- 
cling transition - given the atomic or molecular species 
can be prepared in a cold and slow beam. It may also 
serve as a continuous source of ultra-cold matter that 
may constitute a step towards a continuous wave atom 
laser. 



V. CONCLUSION 



In conclusion, we have developed a model based on 
a the set of coupled rate equations which is capable of 
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